A Multilayered Post-Genome-Wide Association Study Analysis Pipeline Defines Functional Variants and Target Genes for Systemic Lupus Erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWASs) have revealed multiple SLE susceptibility loci and associated single-nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown. METHODS: Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci (cis- and trans-expression quantitative trait loci) with promoter capture high-throughput capture chromatin conformation (PCHi-C), allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, quantitative polymerase chain reaction, and Western blot. RESULTS: Anchoring on 452 index SNPs, we selected 9,931 high linkage disequilibrium (r2 > 0.8) SNPs and defined 182 independent non-human leukocyte antigen (HLA) SLE loci. The 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated reference SNP identifier 57668933 (rs57668933) as a functional variant regulating multiple targets, including SLE-risk gene ELF1 in B cells. CONCLUSION: We demonstrate and validate post-GWAS strategies for using multidimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms to guide experimental characterization.

KLF2 enhancer variant rs4808485 increases lupus risk by modulating inflammasome machinery and cellular homoeostasis.

OBJECTIVE: A recent genome-wide association study linked KLF2 as a novel Asian-specific locus for systemic lupus erythematosus (SLE) susceptibility. However, the underlying causal functional variant(s), cognate target gene(s) and genetic mechanisms associated with SLE risk are unknown. METHODS: We used bioinformatics to prioritise likely functional variants and validated the best candidate with diverse experimental techniques, including genome editing. Gene expression was compared between healthy controls (HCs) and patients with SLE with or without lupus nephritis (LN+, LN-). RESULTS: Through bioinformatics and expression quantitative trait locus analyses, we prioritised rs4808485 in active chromatin, predicted to modulate KLF2 expression. Luciferase reporter assays and chromatin immunoprecipitation-qPCR demonstrated differential allele-specific enhancer activity and binding of active histone marks (H3K27ac, H3K4me3 and H3K4me1), Pol II, CTCF, P300 and the transcription factor PARP1. Chromosome conformation capture-qPCR revealed long-range chromatin interactions between rs4808485 and the KLF2 promoter. These were directly validated by CRISPR-based genetic and epigenetic editing in Jurkat and lymphoblastoid cells. Deleting the rs4808485 enhancer in Jurkat (KO) cells disrupted NLRP3 inflammasome machinery by reducing KLF2 and increasing CASPASE1, IL-1ß and GSDMD levels. Knockout cells also exhibited higher proliferation and cell-cycle progression than wild type. RNA-seq validated interplay between KLF2 and inflammasome machinery in HC, LN+ and LN-. CONCLUSIONS: We demonstrate how rs4808485 modulates the inflammasome and cellular homoeostasis through regulating KLF2 expression. This establishes mechanistic connections between rs4808485 and SLE susceptibility.

Angiogenesis related genes in Takayasu Arteritis (TAK): robust association with Tag SNPs of IL-18 and FGF-2 in a South Asian Cohort.

We performed genetic association study for genes encoding angiogenic and angiostatic proteins in patients with Takayasu arteritis (TAK). A total of 96 SNPs involving 60 genes were studied. Genotyping was performed in Fluidigm 96.96 Dynamic Array chip. All statistical analysis for SNP evaluation was performed using PLINK software. Initial analyses revealed five SNPs from three genes [IL-18 (encodes Interleukin-18), FGF2 (encodes Fibroblast Growth Factor-2), and ANGPT1 (encodes Angiopoietin-1)] as significantly different between controls and cases (uncorrected p < 0.05). After permutation-based analysis, two tag SNPs on the promoter region of IL-18 (rs187238 and rs1946518) and one 3'UTR tag SNP (rs1476217) of FGF2 were significantly associated with susceptibility to TAK, with p and OR (95% CI) of 0.0006 and 1.64 (1.25-2.17), 0.03 and 1.28 (1.02-1.64) & 0.016 and 1.33 (1.05-1.67), respectively; while, the two tag SNPs of ANGPT1 gene (rs6469101 and rs16875900) showed a trend (p = 0.055 & p = 0.051, respectively after permutation based correction). There is robust linkage disequilibrium between the two tag SNPs of IL-18 gene as validated by 1000 genome data of South Asian population; the eQTL effects of these tag SNPs of IL-18 and FGF2 genes on adjacent genes further suggest that these tag SNPs act as genetic risks for development of TAK in South Asians, with possible functional implications towards future biomarker development. Genotype phenotype study by genetic model-based analysis also revealed associations between genotype subsets and clinical features like fever, visual loss, left subclavian and coronary artery involvement in our TAK patients.

A Non-Coding Variant in SLC15A4 Modulates Enhancer Activity and Lysosomal Deacidification Linked to Lupus Susceptibility.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a strong genetic basis. Despite the identification of several single nucleotide polymorphisms (SNPs) near the SLC15A4 gene that are significantly associated with SLE across multiple populations, specific causal SNP(s) and molecular mechanisms responsible for disease susceptibility are unknown. To address this gap, we employed bioinformatics, expression quantitative trait loci (eQTLs), and 3D chromatin interaction analysis to nominate a likely functional variant, rs35907548, in an active intronic enhancer of SLC15A4 . Through luciferase reporter assays followed by chromatin immunoprecipitation (ChIP)-qPCR, we observed significant allele-specific enhancer effects of rs35907548 in diverse cell lines. The rs35907548 risk allele T is associated with increased regulatory activity and target gene expression, as shown by eQTLs and chromosome conformation capture (3C)-qPCR. The latter revealed long-range chromatin interactions between the rs35907548 enhancer and the promoters of SLC15A4, GLTLD1 , and an uncharacterized lncRNA. The enhancer-promoter interactions and expression effects were validated by CRISPR/Cas9 knock-out (KO) of the locus in HL60 promyeloblast cells. KO cells also displayed dramatically dysregulated endolysosomal pH regulation. Together, our data show that the rs35907548 risk allele affects multiple aspects of cellular physiology and may directly contribute to SLE.

Placental microRNA methylome signatures may serve as biomarkers and therapeutic targets for prenatally opioid-exposed infants with neonatal opioid withdrawal syndrome.

Introduction: The neonate exposed to opioids in utero faces a constellation of withdrawal symptoms postpartum commonly called neonatal opioid withdrawal syndrome (NOWS). The incidence of NOWS has increased in recent years due to the opioid epidemic. MicroRNAs (miRNAs) are small non-coding RNA molecules that play a crucial role in gene regulation. Epigenetic variations in microRNAs (miRNAs) and their impact on addiction-related processes is a rapidly evolving area of research. Methods: The Illumina Infinium Methylation EPIC BeadChip was used to analyze DNA methylation levels of miRNA-encoding genes in 96 human placental tissues to identify miRNA gene methylation profiles as-sociated with NOWS: 32 from mothers whose prenatally opioid-exposed infants required pharmacologic management for NOWS, 32 from mothers whose prenatally opioid-exposed infants did not require treat-ment for NOWS, and 32 unexposed controls. Results: The study identified 46 significantly differentially methylated (FDR p-value ≤ 0.05) CpGs associated with 47 unique miRNAs, with a receiver operating characteristic (ROC) area under the curve (AUC) ≥0.75 including 28 hypomethylated and 18 hypermethylated CpGs as potentially associated with NOWS. These dysregulated microRNA methylation patterns may be a contributing factor to NOWS pathogenesis. Conclusion: This is the first study to analyze miRNA methylation profiles in NOWS infants and illustrates the unique role miRNAs might have in diagnosing and treating the disease. Furthermore, these data may provide a step toward feasible precision medicine for NOWS babies as well.

A multilayered post-GWAS analysis pipeline defines functional variants and target genes for systemic lupus erythematosus (SLE).

Objectives: Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWAS) have revealed multiple SLE susceptibility loci and associated single nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown. Methods: Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci ( cis - and trans -eQTLs) with promoter-capture Hi-C, allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, qPCR, and Western blot. Results: Anchoring on 452 index SNPs, we selected 9,931 high-linkage disequilibrium (r 2 >0.8) SNPs and defined 182 independent non-HLA SLE loci. 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated rs57668933 as a functional variant regulating multiple targets, including SLE risk gene ELF1 , in B-cells. Conclusion: We demonstrate and validate post-GWAS strategies for utilizing multi-dimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms, to guide experimental characterization.

Genome-Wide Association Study in Acute Tubulointerstitial Nephritis.

SIGNIFICANCE STATEMENT: Polymorphisms of HLA genes may confer susceptibility to acute tubulointerstitial nephritis (ATIN), but small sample sizes and candidate gene design have hindered their investigation. The first genome-wide association study of ATIN identified two significant loci, risk haplotype DRB1*14-DQA1*0101-DQB1*0503 (DR14 serotype) and protective haplotype DRB1*1501-DQA1*0102-DQB1*0602 (DR15 serotype), with amino acid position 60 in the peptide-binding groove P10 of HLA-DR ß 1 key. Risk alleles were shared among different causes of ATIN and HLA genotypes associated with kidney injury and immune therapy response. HLA alleles showed the strongest association. The findings suggest that a genetically conferred risk of immune dysregulation is part of the pathogenesis of ATIN. BACKGROUND: Acute tubulointerstitial nephritis (ATIN) is a rare immune-related disease, accounting for approximately 10% of patients with unexplained AKI. Previous elucidation of the relationship between genetic factors that contribute to its pathogenesis was hampered because of small sample sizes and candidate gene design. METHODS: We undertook the first two-stage genome-wide association study and meta-analysis involving 544 kidney biopsy-defined patients with ATIN and 2346 controls of Chinese ancestry. We conducted statistical fine-mapping analysis, provided functional annotations of significant variants, estimated single nucleotide polymorphism (SNP)-based heritability, and checked genotype and subphenotype correlations. RESULTS: Two genome-wide significant loci, rs35087390 of HLA-DQA1 ( P =3.01×10 -39 ) on 6p21.32 and rs2417771 of PLEKHA5 on 12p12.3 ( P =2.14×10 -8 ), emerged from the analysis. HLA imputation using two reference panels suggested that HLA-DRB1*14 mainly drives the HLA risk association . HLA-DRB1 residue 60 belonging to pocket P10 was the key amino acid position. The SNP-based heritability estimates with and without the HLA locus were 20.43% and 10.35%, respectively. Different clinical subphenotypes (drug-related or tubulointerstitial nephritis and uveitis syndrome) seemed to share the same risk alleles. However, the HLA risk genotype was associated with disease severity and response rate to immunosuppressive therapy. CONCLUSIONS: We identified two candidate genome regions associated with susceptibility to ATIN. The findings suggest that a genetically conferred risk of immune dysregulation is involved in the pathogenesis of ATIN.

Placental cytochrome P450 methylomes in infants exposed to prenatal opioids: exploring the effects of neonatal opioid withdrawal syndrome on health horizons.

Background: Neonatal opioid withdrawal syndrome (NOWS), arises due to increased opioid use during pregnancy. Cytochrome P450 (CYP) enzymes play a pivotal role in metabolizing a wide range of substances in the human body, including opioids, other drugs, toxins, and endogenous compounds. The association between CYP gene methylation and opioid effects is unexplored and it could offer promising insights. Objective: To investigate the impact of prenatal opioid exposure on disrupted CYPs in infants and their anticipated long-term clinical implications. Study Design: DNA methylation levels of CYP genes were analyzed in a cohort of 96 placental tissues using Illumina Infinium MethylationEPIC (850 k) BeadChips. This involved three groups of placental tissues: 32 from mothers with infants exposed to opioids prenatally requiring pharmacologic treatment for NOWS, 32 from mothers with prenatally opioid-exposed infants not needing NOWS treatment, and 32 from unexposed control mothers. Results: The study identified 20 significantly differentially methylated CpG sites associated with 17 distinct CYP genes, with 14 CpGs showing reduced methylation across 14 genes (CYP19A1, CYP1A2, CYP4V2, CYP1B1, CYP24A1, CYP26B1, CYP26C1, CYP2C18, CYP2C9, CYP2U1, CYP39A1, CYP2R1, CYP4Z1, CYP2D7P1 and), while 8 exhibited hypermethylation (CYP51A1, CYP26B1, CYP2R1, CYP2U1, CYP4X1, CYP1A2, CYP2W1, and CYP4V2). Genes such as CYP1A2, CYP26B1, CYP2R1, CYP2U1, and CYP4V2 exhibited both increased and decreased methylation. These genes are crucial for metabolizing eicosanoids, fatty acids, drugs, and diverse substances. Conclusion: The study identified profound methylation changes in multiple CYP genes in the placental tissues relevant to NOWS. This suggests that disruption of DNA methylation patterns in CYP transcripts might play a role in NOWS and may serve as valuable biomarkers, suggesting a future pathway for personalized treatment. Further research is needed to confirm these findings and explore their potential for diagnosis and treatment.

A Non-Coding Variant in SLC15A4 Modulates Enhancer Activity and Lysosomal Deacidification Linked to Lupus Susceptibility.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a strong genetic basis. Despite the identification of several single nucleotide polymorphisms (SNPs) near the SLC15A4 gene that are significantly associated with SLE across multiple populations, specific causal SNP(s) and molecular mechanisms responsible for disease susceptibility are unknown. To address this gap, we employed bioinformatics, expression quantitative trait loci (eQTLs), and 3D chromatin interaction analysis to nominate a likely functional variant, rs35907548, in an active intronic enhancer of SLC15A4. Through luciferase reporter assays followed by chromatin immunoprecipitation (ChIP)-qPCR, we observed significant allele-specific enhancer effects of rs35907548 in diverse cell lines. The rs35907548 risk allele T is associated with increased regulatory activity and target gene expression, as shown by eQTLs and chromosome conformation capture (3C)-qPCR. The latter revealed long-range chromatin interactions between the rs35907548 enhancer and the promoters of SLC15A4, GLTLD1, and an uncharacterized lncRNA. The enhancer-promoter interactions and expression effects were validated by CRISPR/Cas9 knock-out (KO) of the locus in HL60 promyeloblast cells. KO cells also displayed dramatically dysregulated endolysosomal pH regulation. Together, our data show that the rs35907548 risk allele affects multiple aspects of cellular physiology and may directly contribute to SLE.

Genome-wide association study for systemic lupus erythematosus in an egyptian population.

Systemic lupus erythematosus (SLE) susceptibility has a strong genetic component. Genome-wide association studies (GWAS) across trans-ancestral populations show both common and distinct genetic variants of susceptibility across European and Asian ancestries, while many other ethnic populations remain underexplored. We conducted the first SLE GWAS on Egyptians-an admixed North African/Middle Eastern population-using 537 patients and 883 controls. To identify novel susceptibility loci and replicate previously known loci, we performed imputation-based association analysis with 6,382,276 SNPs while accounting for individual admixture. We validated the association analysis using adaptive permutation tests (n = 109). We identified a novel genome-wide significant locus near IRS1/miR-5702 (Pcorrected = 1.98 × 10-8) and eight novel suggestive loci (Pcorrected < 1.0 × 10-5). We also replicated (Pperm < 0.01) 97 previously known loci with at least one associated nearby SNP, with ITGAM, DEF6-PPARD and IRF5 the top three replicated loci. SNPs correlated (r 2 > 0.8) with lead SNPs from four suggestive loci (ARMC9, DIAPH3, IFLDT1, and ENTPD3) were associated with differential gene expression (3.5 × 10-95 < p < 1.0 × 10-2) across diverse tissues. These loci are involved in cellular proliferation and invasion-pathways prominent in lupus and nephritis. Our study highlights the utility of GWAS in an admixed Egyptian population for delineating new genetic associations and for understanding SLE pathogenesis.

Discovery and Functional Characterization of Two Regulatory Variants Underlying Lupus Susceptibility at 2p13.1.

Genome-wide association studies have identified 2p13.1 as a prominent susceptibility locus for systemic lupus erythematosus (SLE)­a complex, multisystem autoimmune disease. However, the identity of underlying causal variant (s) and molecular mechanisms for increasing disease susceptibility are poorly understood. Using meta-analysis (cases = 10,252, controls = 21,604) followed by conditional analysis, bioinformatic annotation, and eQTL and 3D-chromatin interaction analyses, we computationally prioritized potential functional variants and subsequently experimentally validated their effects. Ethnicity-specific meta-analysis revealed striking allele frequency differences between Asian and European ancestries, but with similar odds ratios. We identified 20 genome-wide significant (p < 5 × 10−8) variants, and conditional analysis pinpointed two potential functional variants, rs6705628 and rs2272165, likely to explain the association. The two SNPs are near DGUOK, mitochondrial deoxyguanosine kinase, and its associated antisense RNA DGUOK-AS1. Using luciferase reporter gene assays, we found significant cell type- and allele-specific promoter activity at rs6705628 and enhancer activity at rs2272165. This is supported by ChIP-qPCR showing allele-specific binding with three histone marks (H3K27ac, H3K4me3, and H3K4me1), RNA polymerase II (Pol II), transcriptional coactivator p300, CCCTC-binding factor (CTCF), and transcription factor ARID3A. Transcriptome data across 28 immune cell types from Asians showed both SNPs are cell-type-specific but only in B-cells. Splicing QTLs showed strong regulation of DGUOK-AS1. Genotype-specific DGOUK protein levels are supported by Western blots. Promoter capture Hi-C data revealed long-range chromatin interactions between rs2272165 and several nearby promoters, including DGUOK. Taken together, we provide mechanistic insights into how two noncoding variants underlie SLE risk at the 2p13.1 locus.

Genetics Matters: Voyaging from the Past into the Future of Humanity and Sustainability.

The understanding of how genetic information may be inherited through generations was established by Gregor Mendel in the 1860s when he developed the fundamental principles of inheritance. The science of genetics, however, began to flourish only during the mid-1940s when DNA was identified as the carrier of genetic information. The world has since then witnessed rapid development of genetic technologies, with the latest being genome-editing tools, which have revolutionized fields from medicine to agriculture. This review walks through the historical timeline of genetics research and deliberates how this discipline might furnish a sustainable future for humanity.

Maternal opioid use disorder: Placental transcriptome analysis for neonatal opioid withdrawal syndrome.

Excessive prenatal opioid exposure may lead to the development of Neonatal Opioid Withdrawal Syndrome (NOWS). RNA-seq was done on 64 formalin-fixed paraffin-embedded placental tissue samples from 32 mothers with opioid use disorder, with newborns with NOWS that required treatment, and 32 prenatally unexposed controls. We identified 93 differentially expressed genes in the placentas of infants with NOWS compared to unexposed controls. There were 4 up- and 89 downregulated genes. Among these, 7 genes CYP1A1, APOB, RPH3A, NRXN1, LINC01206, AL157396.1, UNC80 achieved an FDR p-value of <0.01. The remaining 87 genes were significant with FDR p-value <0.05. The 4 upregulated, CYP1A1, FP671120.3, RAD1, RN7SL856P, and the 10 most significantly downregulated genes were RNA5SP364, GRIN2A, UNC5D, DMBT1P1, MIR3976HG, LINC02199, LINC02822, PANTR1, AC012178.1, CTNNA2. Ingenuity Pathway Analysis identified the 7 most likely to play an important role in the etiology of NOWS. Our study expands insights into the genetic mechanisms of NOWS development.

Biallelic variants in KARS1 are associated with neurodevelopmental disorders and hearing loss recapitulated by the knockout zebrafish.

PURPOSE: Pathogenic variants in Lysyl-tRNA synthetase 1 (KARS1) have increasingly been recognized as a cause of early-onset complex neurological phenotypes. To advance the timely diagnosis of KARS1-related disorders, we sought to delineate its phenotype and generate a disease model to understand its function in vivo. METHODS: Through international collaboration, we identified 22 affected individuals from 16 unrelated families harboring biallelic likely pathogenic or pathogenic in KARS1 variants. Sequencing approaches ranged from disease-specific panels to genome sequencing. We generated loss-of-function alleles in zebrafish. RESULTS: We identify ten new and four known biallelic missense variants in KARS1 presenting with a moderate-to-severe developmental delay, progressive neurological and neurosensory abnormalities, and variable white matter involvement. We describe novel KARS1-associated signs such as autism, hyperactive behavior, pontine hypoplasia, and cerebellar atrophy with prevalent vermian involvement. Loss of kars1 leads to upregulation of p53, tissue-specific apoptosis, and downregulation of neurodevelopmental related genes, recapitulating key tissue-specific disease phenotypes of patients. Inhibition of p53 rescued several defects of kars1-/- knockouts. CONCLUSION: Our work delineates the clinical spectrum associated with KARS1 defects and provides a novel animal model for KARS1-related human diseases revealing p53 signaling components as potential therapeutic targets.

Lupus Susceptibility Region Containing CDKN1B rs34330 Mechanistically Influences Expression and Function of Multiple Target Genes, Also Linked to Proliferation and Apoptosis.

OBJECTIVE: In a recent genome-wide association study, a significant genetic association between rs34330 of CDKN1B and risk of systemic lupus erythematosus (SLE) in Han Chinese was identified. This study was undertaken to validate the reported association and elucidate the biochemical mechanisms underlying the effect of the variant. METHODS: We performed an allelic association analysis in patients with SLE, followed by a meta-analysis assessing genome-wide association data across 11 independent cohorts (n = 28,872). In silico bioinformatics analysis and experimental validation in SLE-relevant cell lines were applied to determine the functional consequences of rs34330. RESULTS: We replicated a genetic association between SLE and rs34330 (meta-analysis P = 5.29 × 10-22 , odds ratio 0.84 [95% confidence interval 0.81-0.87]). Follow-up bioinformatics and expression quantitative trait locus analysis suggested that rs34330 is located in active chromatin and potentially regulates several target genes. Using luciferase and chromatin immunoprecipitation-real-time quantitative polymerase chain reaction, we demonstrated substantial allele-specific promoter and enhancer activity, and allele-specific binding of 3 histone marks (H3K27ac, H3K4me3, and H3K4me1), RNA polymerase II (Pol II), CCCTC-binding factor, and a critical immune transcription factor (interferon regulatory factor 1 [IRF-1]). Chromosome conformation capture revealed long-range chromatin interactions between rs34330 and the promoters of neighboring genes APOLD1 and DDX47, and effects on CDKN1B and the other target genes were directly validated by clustered regularly interspaced short palindromic repeat (CRISPR)-based genome editing. Finally, CRISPR/dead CRISPR-associated protein 9-based epigenetic activation/silencing confirmed these results. Gene-edited cell lines also showed higher levels of proliferation and apoptosis. CONCLUSION: Collectively, these findings suggest a mechanism whereby the rs34330 risk allele (C) influences the presence of histone marks, RNA Pol II, and IRF-1 transcription factor to regulate expression of several target genes linked to proliferation and apoptosis. This process could potentially underlie the association of rs34330 with SLE.

Deep sequencing reveals a DAP1 regulatory haplotype that potentiates autoimmunity in systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a clinically heterogeneous autoimmune disease characterized by the development of anti-nuclear antibodies. Susceptibility to SLE is multifactorial, with a combination of genetic and environmental risk factors contributing to disease development. Like other polygenic diseases, a significant proportion of estimated SLE heritability is not accounted for by common disease alleles analyzed by SNP array-based GWASs. Death-associated protein 1 (DAP1) was implicated as a candidate gene in a previous familial linkage study of SLE and rheumatoid arthritis, but the association has not been explored further. RESULTS: We perform deep sequencing across the DAP1 genomic segment in 2032 SLE patients, and healthy controls, and discover a low-frequency functional haplotype strongly associated with SLE risk in multiple ethnicities. We find multiple cis-eQTLs embedded in a risk haplotype that progressively downregulates DAP1 transcription in immune cells. Decreased DAP1 transcription results in reduced DAP1 protein in peripheral blood mononuclear cells, monocytes, and lymphoblastoid cell lines, leading to enhanced autophagic flux in immune cells expressing the DAP1 risk haplotype. Patients with DAP1 risk allele exhibit significantly higher autoantibody titers and altered expression of the immune system, autophagy, and apoptosis pathway transcripts, indicating that the DAP1 risk allele mediates enhanced autophagy, leading to the survival of autoreactive lymphocytes and increased autoantibody. CONCLUSIONS: We demonstrate how targeted sequencing captures low-frequency functional risk alleles that are missed by SNP array-based studies. SLE patients with the DAP1 genotype have distinct autoantibody and transcription profiles, supporting the dissection of SLE heterogeneity by genetic analysis.

ITGAM is a risk factor to systemic lupus erythematosus and possibly a protection factor to rheumatoid arthritis in patients from Mexico.

INTRODUCTION: ITGAM has consistently been associated with susceptibility to systemic lupus erythematosus (SLE) in many ethnically diverse populations. However, in populations with higher Amerindian ancestry (like Yucatan) or highly admixed population (like Mexican), ITGAM has seldom been evaluated (except few studies where patients with childhood-onset SLE were included). In addition, ITGAM has seldom been evaluated in patients with rheumatoid arthritis (RA). Here, we evaluated whether four single nucleotide polymorphisms (SNPs), located within ITGAM, were associated with SLE and RA susceptibility in patients from Mexico. METHODS: Our study consisted of 1,462 individuals, which included 363 patients with SLE (292 from Central Mexico and 71 from Yucatan), and 621 healthy controls (504 from Central Mexico and 117 from Yucatan). Our study also included 478 patients with RA from Central Mexico. TaqMan assays were used to obtain the genotypes of the four SNPs: rs34572943 (G/A), rs1143679 (G/A), rs9888739 (C/T), and rs1143683 (C/T). We also verified the genotypes by Sanger sequencing. Fisher's exact test and permutation test were employed to evaluate the distribution of genotype, allele, and haplotype between patients and controls. RESULTS: Our data show that all four ITGAM SNPs are significantly associated with susceptibility to SLE using both genotypic and allelic association tests (corrected for multiple testing, but not for population stratification). A second study carried out in patients from Yucatan, a southeastern part of Mexico (with a high Amerindian ancestry), also replicated SLE association with all four SNPs, including the functional SNP, rs1143679 (OR = 24.6 and p = 9.3X10-6). On the other hand, none of the four SNPs are significant in RA after multiple testing. Interestingly, the GACC haplotype, which carries the ITGAM rs1143679 (A) minor allele, showed an association with protection against RA (OR = 0.14 and p = 3.0x10-4). CONCLUSION: Our data displayed that ITGAM is a risk factor to SLE in individuals of Mexican population. Concurrently, a risk haplotype in ITGAM confers protection against RA.

Mechanistic Characterization of RASGRP1 Variants Identifies an hnRNP-K-Regulated Transcriptional Enhancer Contributing to SLE Susceptibility.

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component. We recently identified a novel SLE susceptibility locus near RASGRP1, which governs the ERK/MAPK kinase cascade and B-/T-cell differentiation and development. However, precise causal RASGRP1 functional variant(s) and their mechanisms of action in SLE pathogenesis remain undefined. Our goal was to fine-map this locus, prioritize genetic variants likely to be functional, experimentally validate their biochemical mechanisms, and determine the contribution of these SNPs to SLE risk. We performed a meta-analysis across six Asian and European cohorts (9,529 cases; 22,462 controls), followed by in silico bioinformatic and epigenetic analyses to prioritize potentially functional SNPs. We experimentally validated the functional significance and mechanism of action of three SNPs in cultured T-cells. Meta-analysis identified 18 genome-wide significant (p < 5 × 10-8) SNPs, mostly concentrated in two haplotype blocks, one intronic and the other intergenic. Epigenetic fine-mapping, allelic, eQTL, and imbalance analyses predicted three transcriptional regulatory regions with four SNPs (rs7170151, rs11631591-rs7173565, and rs9920715) prioritized for functional validation. Luciferase reporter assays indicated significant allele-specific enhancer activity for intronic rs7170151 and rs11631591-rs7173565 in T-lymphoid (Jurkat) cells, but not in HEK293 cells. Following up with EMSA, mass spectrometry, and ChIP-qPCR, we detected allele-dependent interactions between heterogeneous nuclear ribonucleoprotein K (hnRNP-K) and rs11631591. Furthermore, inhibition of hnRNP-K in Jurkat and primary T-cells downregulated RASGRP1 and ERK/MAPK signaling. Comprehensive association, bioinformatics, and epigenetic analyses yielded putative functional variants of RASGRP1, which were experimentally validated. Notably, intronic variant (rs11631591) is located in a cell type-specific enhancer sequence, where its risk allele binds to the hnRNP-K protein and modulates RASGRP1 expression in Jurkat and primary T-cells. As risk allele dosage of rs11631591 correlates with increased RASGRP1 expression and ERK activity, we suggest that this SNP may underlie SLE risk at this locus.

Amino acid signatures of HLA Class-I and II molecules are strongly associated with SLE susceptibility and autoantibody production in Eastern Asians.

Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10-75) and its proxy position 11 (P = 1.1x10-67), followed by HLA-DRB1-37 (P = 4.5x10-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10-8), HLA-DPB1-35 (P = 9.0x10-16), HLA-DQB1-37 (P = 2.7x10-14), and HLA-B-9 (P = 6.5x10-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10-14; DRB1-13, P = 2.9x10-13; DRB1-30, P = 3.9x10-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.

Identification of Proteins Interacting with Single Nucleotide Polymorphisms (SNPs) by DNA Pull-Down Assay.

Single nucleotide polymorphisms (SNPs) are the most abundantly found common form of DNA variation in the human genome. Many genetic association studies have proved that some of these SNPs are involved in regulating several cellular/physiological processes ranging from gene regulation to disease development. Analysis of the protein complex that binds to these SNPs is a crucial step in studying the mechanisms by which gene expressions are regulated in cis- or trans-acting manner. Commonly used techniques to determine DNA-protein interaction, such as electrophoretic mobility shift assay (EMSA), have limited value for simultaneously analyzing a large number of proteins in the complex. Furthermore, this assay is tedious and time-consuming and often requires radiolabeled probe as well as extensive optimization. Here, we describe a pull-down assay before performing the EMSA, which helps in the detection of differentially-bound protein(s) in an allele-specific manner. The assay is easy to perform and does not require radiolabeling of DNA probes. Biotinylated DNA probe bound to streptavidin beads can be complexed with protein(s) from cell nuclear lysate, nonspecific proteins were washed out, and only protein(s) having high affinity to SNP-specific DNA were detected on SDS-PAGE and identified by mass spectrometry.